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Abstract 

We analyse the possible existence of non-perturbative contributions in heavy 
QQ systems (Q and Q need not have the same flavour) which cannot be ex- 
pressed in terms of local condensates. Starting from QCD, with well defined 
approximations and splitting properly the fields into large and small momen- 
tum components, we derive an effective lagrangian where hard gluons (in the 
non-relativistic aproximation) have been integrated out. The large momentum 
contributions (which are dominant) are calculated using Coulomb type states. 
Besides the usual condensate corrections, we see the possibility of new non- 
perturbative contributions. We parametrize them in terms of two low momentum 
correlators with Coulomb bound state energy insertions En- We realize that the 
Heavy Quark Effective lagrangian can be used in these correlators. We calculate 
the corrections that they give rise to in the decay constant, the bound state en- 
ergy and the matrix elements of bilinear currents at zero recoil. We study the 
cut-off dependence of the new contributions and we see that it matches perfectly 
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with that of the large momentum contributions. We consider two situations in 
detail: i) En » ^qcd {Mq — oo) and ii) En « ^qcd, and briefly discuss 
the expected size of the new contributions in T , J/^* and B* systems. 
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1 Introduction 



The study of heavy quark bound state systems remains one of the more promising 
topics in order to test both perturbative and non-perturbative aspects of QCD, as it is 
clear from the steady activity in the field [0-0 • These systems can be understood in 
a first approximation as non-relativistic bound states which occur due to a Coulomb 
type interaction predicted by perturbative QCD. In order to improve this basic picture 
one has to deal on one side with perturbative relativistic and radiative corrections, and 
on the other side with non-perturbative corrections (power corrections). 

In this paper we shall only be concerned with non-perturbative corrections. Usu- 
ally, the latter have been parametrized using both the multipole expansion and the 
adiabatic approximation in terms of the gluon condensate 0, §]. Corrections to the 
Coulomb potential due to condensates can also be considered, although these are sub- 
leading ||. We have argued before [§] that new non-perturbative contributions could 
arise which cannot be expressed in terms of local condensates, and hence a convenient 
parametrization for them is required. This kind of nonperturbative contributions has 
been discussed in |]10[ in a different context and, in fact, the various Isgur-Wise func- 
tions extensively used in the Heavy Quark Effective Theory (HQET) may be regarded 
as such [|lT] . 

Let us recall the main idea behind the possibility of new non-perturbative contri- 
butions in heavy quarkoniumQ. When the relative three momentum in the bound state 
is big enough the dominant interaction must be the perturbative Coulomb potential, 
but for small relative three momentum this need not longer be true. Therefore, heavy 
quarks in the latter kinematical situation should better be kept as low energy degrees 
of freedom. It turns out that a convenient parametrization of this kinematical region 



may be given in terms of the HQET for quarks and antiquarks |T2|. 

The HQET for quarks and antiquarks enjoys rather peculiar features, which make 
it quite different from the usual HQET describing either quarks or antiquarks, which 
has been so popular in the study of Qq and Qqq systems in recent years [O (see WA 



for reviews). For instance, it enjoys a symmetry, which is larger than the well-known 
spin and flavor symmetry, that breaks spontanously down to the latter giving rise to 



"'^We use 'heavy quarkonium' to denote a general heavy quark-antiquark bound state. The quark 
and the antiquark need not have the same flavour. 
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quark antiquark states as Goldstone modes [|T2|. Its pecularities concerning radiative 



corrections have recently been illustrated in ||T5|. 

The main aim of this paper is to work out a controlled derivation from QCD of the 
effective lagrangian describing the small relative momentum regime of heavy quarks 
in quarkonium. Whereas the basic ideas above have already been elaborated in 0, 
a complete and systematic derivation is lacking, and hence worth being presented. 
Within this new framework we recalculate the non-perturbative contributions of this 
region to the energy levels, the decay constant and the matrix elements of bilinear 
currents at zero recoil. We find a few corrections to the formulas given in 0. For all 
these observables it is enough to work in the center of mass frame (CM), which we 
shall do in most of the paper. 

In order to deal with heavy quarkonia systems we keep the relevant degrees of 
freedom in the QCD lagrangian. In fact, since virtual heavy quark creation is very 
much suppressed, we could safely start from non-relativistic QCD (NRQCD). The 
derivation of NRQCD from QCD is well understood and a technique to incorporate 
relativistic corrections to it has also been developed |]T6[. First of all, we split the gluon 
field in hard and soft by a three momentum cut-off. From the hard gluon fields we 
only keep the zero component and disregard the spacial components. This is legitimate 
as far as we are not interested in relativistic corrections. We next integrate out the 
zero component of the hard gluon field to obtain the Coulomb potential between heavy 
quark currents. The Coulomb potential has an infrared momentum cut-off since the 
zero component of the soft gluon field has not been integrated out. At this point we have 
an effective lagrangian formally equal to the one used by Voloshin and Leutwyler (VL) 
U^, §], except for the IR cut-off in the Coulomb potential. After introducing CM and 
relative momentum for the bound states we are interested in, we further split the quark 
fields in large and small relative three-momentum regimes^. The resulting lagrangian 
can then be separated in three pieces: L'^ which contains small relative momentum 
quark fields only, + contains large relative momentum quark fields only and 
which contains both small and large relative momentum quark fields. For L^^ we can 
aproximate the lagrangian to the HQET lagrangian, where eventually all its powerful 
symmetries can be used. No Coulomb term remains in this part of the lagrangian. For 

+ we obtain again the VL starting point lagrangian except for two facts: both 

^The large and small relative momentum regions were denoted as off- and on-shell regions respec- 
tively in [§. 
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the Coulomb potential and the Hilbert space are restricted to three-momenta larger 
than a certain cut-off. Keeping the cut-off much higher than Aqcd but much smaller 
than the invers Bohr radius we may safely assume that the multipole expansion holds 
for this part of the lagrangian. If we further assume that the adiabatic approximation 
also holds, we may proceed in total analogy to VL. The hipothesis above on the cut- 
off also allows us to treat L^^ as a perturbation. The various contributions from this 
perturbation to the different observables can be eventually expressed as correlators of 
the HQET. 

We would hke to stress that our formahsm is less restrictive than the one used by 
VL since neither the adiabatic approximation nor the multipole expansion are assumed 
to hold in the small relative momentum region of the heavy quark fields. Indeed we 
may always recover VLs results by putting to zero the cut-off which separates large 
and small relative momentum. 

We distribute the paper as follows. In sect. 2 we derive the effective action for 
the small relative momentum fields. In sect. 3 we calculate the decay constant, the 
bound state mass and the matrix element of any bilinear heavy quark current between 
quarkonia states at zero recoil. The latter is relevant for the study of semileptonic 
decays at zero recoil. In sect. 4 we prove the cut-off independence of our results. In 
section 5 we study the low momentum correlators in two situations: the asymptotic 
limit {Mq — > oo) En » Aqcd, where, using OPE techniques, we see that no new 
corrections arise, and ii) En « Aqcd, where the low momentum contributions are 
evaluated using an effective 'chiral' lagrangian which incorporates the relevant sym- 
metries of the HQET for quarks and antiquarks. Working in this way we find new 
non-perturbative contributions which are parametrized by a single non-perturbative 
constant. We also give preliminary estimations of their size. Section 6 is devoted to 
the conclusions. 

2 Effective action 

In this section we derive the effective lagrangian for heavy quarks and antiquarks 
in the small relative momentum regime from QCD. 
The QCD lagrangian reads 

^--\f^ + J2 Qa{ip- ma)Qa (2.1) 
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where 

D^ = d^~igV^, V = VT'- (2.2) 

f;. = d,v: - d.v; + gr'vx (2-3) 

We split the gluon field V in large A and small B momentum modes V{x) = 
A{x) + B{x). Next we exactly integrate Ao and neglect Ai. The latter would give 
rise to relativistic corrections. Consistently, at the same point we perform a Foldy- 
Wouthuysen transformation and keep terms up to 1/m. We obtain 

L = -1 1 £xFl + Y.j [Qaill^'D^ - ma)Qa + Qa^Qa + Q-^^^Q- 

1 2 / 1 \ 

- yE/^'^/ d'yQal'T^Qaix) {x,y)Qa'l'T'Qa'{y) (2.4) 

a aa' \ ^ B ' 

which is manifestly gauge invariant ^ . Although, in principle, we could attempt to 
carry out an explicitely gauge invariant calculation, in practise, it is most convenient 
to choose a slightly modified Schwinger gauge for the small momentum gluons 

{z )5(z) = (2.5) 

ma + vrta' 

In this gauge B in the kinetic and Coulomb terms gives rise to subdominant contribu- 
tions when the multipole expansion is carried out, which greatly simplifies the calcula- 
tion. In particular, recall that the propagator in the Coulomb term always carries large 
momentum (we have not integrated out the small momentum Vq which is kept in B^. 
Hence the multipole expansion is always legitimated in the Coulomb term. This allows 
to drop B in the Coulomb term straight away. As long as we are interested in quark- 
antiquark bound states only, we may also safely neglect the four-fermion interaction 
terms involving only quarks or only antiquarks. We next rearrange the quark- ant iquark 
interaction term in a convenient way in order to describe the bound state dynamics . 
Finally, the effective lagrangian reads 

L = -1 1 d'xFl +Y,j£x |^Qa(^7°< - ma)Qa + Qa^Qa 

X Q^{-maV + q,t)T'^T,Qa>{'maiV + q,t) Qa>{ma'V + ^ ,t)T^T,Qai-maV + g' ,t) 



^Similar approaches can be found in the literature [Q. 
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where A = 0,r denotes colour (0 singlet and r octet, r = 1...8), \q — q^l > n , n being 
the cut-off which separates small and large momenta, and 

1±70 

rriaa' =ma + ma' , = «75p_ , i^P-, P± ■= — n — ' (2-7) 



While, the potential reads 



V%P) = -^, V^{p) = -^. (2.8) 

where Cp = {N^ — l)/2A'c and |p| > must be understood due to the cut-off coming 
from soft gluons0. 

Written in this way, we can understand the four-fermion Coulomb interaction term 
as one which creates a quark-antiquark state with central velocity v and relative mo- 
mentum (f and annihilates a quark-antiquark state with the same center of masses 
velocity v and relative momentum q. Obviously v is a. conserved quantity in this 
non-relativistic approximation. We consider the spin breaking term as subleading and 
we will neglect it in the following. Therefore, spin symmetry for both low and high 
momentum is implicit in the rest of the paper. 

If we stopped at this point we would obtain the standard VL results. However, we 
would like to go beyond and look for new non-perturbative contributions. We observe 
that quarks with small relative three momentum only feel the Coulomb interaction 
of quarks with large relative momentum. This suggests to perform a splitting of the 
physical quark and antiquark fields into small and large relative momentum in the 
bound state. The physical picture behind is that if the relative three momentum in 
the bound state is big enough we can understand it as a perturbative Coulomb type 
bound state. But for small relative three momentum that is no longer true. For that 
momentum regime the quark and antiquark fields should be kept as low momentum 
degrees of freedom. That is, in fact, the main idea of the paper. Therefore, let us write 
down the currents related to the physical quark-antiquark states in momentum space 

I (0^'"""'"''''7 (03§a'(-"^a'^- g,t)T^rQaKt/- q,t) . (2.9) 



^Several aspects related to this cut-ofF dependence have been studied in pi 
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The matrix T should be such that it projects over quark-antiquark states according to 
our non-relativistic picture. Notice that the time dependence is kept exphcit. Further- 
more, we spht the relative three momentum with the same cut-off as above. Thus, 
( ^ ) reads 

Jp {X) — Y \X) + p y^x) 



+ J (2^^'"""'" 7;. j^^QA-^a'V-q,t)T''TQ,{m,v- q,t) , (2.10) 

where Qaif^aV + k) =: h]!^{k). After that we may divide the lagrangian in three pieces. 

1 = 1^ + 1^^ + 1' (2.11) 
is the piece of the effective lagrangian containing large momenta only. It reads 



^t^ = Y.j ^Qaih^do - ma)Qa + Q<^^Qo^ 



^ aa'A 



r d?v 


f <fq' f 


1 (27r)3 


i (27r)3 i 



(2.12) 



X 



[Q^{-maV + q,t)T'^TsQa'{'ma'V + q,t) Q^'irria^v + q' ,t)T'^T sQa{-maV + q',t) 



where \q — (f \ > fi. 

In fact it is nothing but the standard Coulomb lagrangian, except for the cut-offs, 
is the piece of the effective lagrangian containing small momenta only. It reads 

= ~\fI + E \h::{il'D^ - m.)K + K^k\^ ■ (2.13) 

Notice that ( |2.13| ) does not have the four-fermion Coulomb term. It contains the whole 
soft gluon lagrangian as well as the heavy quark and antiquark fields with small three 
relative momentum. All the fields in (|2.13|) are in the non-perturbative regime of QCD. 
Notice that if we drop the term in l/m^ and make K"^ g-n^^iaxo^?; ( |2.13| ) becomes 
the HQET lagrangian in the rest frame. Although the 1 / ma term is naively subleading 
for small relative momentum, it plays a crucial role in certain circumstances, as we 
shall see in Section 4. Nonetheless let us advance that for the correlators we will be 
interested in one can safely neglect it and work with the HQET lagrangian. 
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U mixes small and large momenta 

L' = Li + Li^ (2.14) 



The first term reads 



Ci{x)=gQal'BirQa{x) (2.15) 

which gives the leading contribution to the multipole expansion. We will not discuss 
these contributions ( p.l5| ) here since they have been extensively studied in the literature 



1 1^, . Let us focus on the second term. It reads 



L-^-\ E r^.Nl ^j^^l^f ^V^^,^ - ,T ,2.16) 



aa'A 



Ji\{q,t)T^T,hl,{q,t)] [^,,(m,,t; + g', t)r^r,g,(-m,t; + g', t)] + (h.c.) . 

In this expression the Coulomb potential is the only piece which mixes small and large 
relative momentum. We can perform a derivative expansion since q and q' belong to 
different momentum regimes {q ~ Aqcd << <?' ~ ~) and keep only the leading term 
(further orders would give subleading corrections). It turns out that the small relative 
momentum term decouples from the Coulomb potential and can be written like a local 
current. Finally, we obtain 

X J d^xe-'^^^"'-^Jl^^;'{x,t) [Q,,{ma'V + q,t)T^rsQa{-maV + q,t)]+{h.c.). 

The formalism developed in was not powerful enough as to uncover the in- 
teraction lagrangian (2.17). This interaction lagrangian is indeed the responsible for 
the differences between the results presented there and the ones obtained in the next 
section. 

If we assume that small momentum terms are small in comparison with the large 
momentum terms we can treat the interaction lagrangian (2.17) as a perturbation. 
This is so for the lower energy levels of heavy quark bound states. In the next sections 
we focus on the nonperturbative contributions coming from (2.13) and (2.17). 
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3 Physical Observables 

In this section we work out the non-perturbative corrections from the small relative 
momentum region to the decay constant, the bound state mass and the matrix elements 
of bilinear currents at zero recoil. We take the bound state velocity small or zero. 

Consider first the eigenvalues and eigenstates of H^, the hamiltonian associated to 
L^. They read 

\{ah,n,s,A)-v) = -| j ^^.^Xn^t, i^)u-{p,)V ,v%p,^^^^^^ , 

(3.1) 

EiM (3.2) 

where 

— * — * 

Pi=maV + k, p2 = mbV-k, (3.3) 

Eab,ny ^ab,ni^) ^ab,ni^) ^he energy, the coordinate space wave function and 
the momentum space wave function of a Coulomb-type state with quantum number 
n = {n,l,m). v is the bound state 3- vector velocity, a and b are flavour indices 
and s denotes spin, b'^ and S are creation operators of particles and anti-particles 
respectively. 

{bt^iPi)^ bi,M)} = {2irfS^%pS,,,,6'{p^ - p,) (3.4) 

<..(P2)} = {2Txf5^'5.p5,,,,5\p, - p2) (3.5) 

{<n(Pi)'<..(P2)} = (3.6) 

u°'{pi) and v^{p2) are spinors normalized in such a way that in the large m limit 
the following holds 

u'^{pi)u^{pi) = p+ , E = -P- ■ (3-7) 

a a 

( |3.1D has the non-relativistic normalization 

{{ab,n,s,A);v\{a'b',n',s',A');v') = {27rfS^^\mabiv - v'))5n,n'5s,s'S(ab),{a'b')SA,A' , (3.8) 



where we have used 

tr {p+rp_r') = -25s,s' 



(3.9) 
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From ( |3.8D it follows the wave function normalization 



where there is no sum over A. The wave function and the energy fulfil the equation 



marrib 

p > /i , fXab = ■ ■ 

rua + rub 

From (|2.^ ) it trivially follows that the eight components of the octet wave function 
fulfil the same equation and hence they are the same. Notice that the wave function 
normalization and the differential equation above are fi dependent. Furthermore, the 
wave function is not defined over all values of p. We will work this out in detail in 
sect. 4. In order to simplify the notation we will not displayed the cut-off dependence 
explicitely in the rest of the section, but it must be understood throughout. 

For H'^, the hamiltonian associated to L^, we denote the eigenstates and eigenvalues 

by 

\{ab,g,s);v), Eg (3.12) 

where g labels the low momentum state. We cannot give explicit expressions since their 
dynamics is governed by low momentum. (|3.12| ) has the non-relativistic normalization 

{{ab,g,s);v\{a'b',g',s');v') = i27rYS^^\mab{v-v))Ss,s'S(ab),ia'b')Sg,g'- (3.13) 

Of course the states ( |3Tl| ) and ( |3.12| ) are orthogonal since they belong to different 
momentum regimes. 

Our Hilbert space is (before switching on Ljf-) {N}={(n,A), g} and the identity 
reads in this base 

r d^P 

1^|0)(0| + P' + V=|0)(0|+ Yl J -7^\{ab,N,s);v){{ab,N,s);v\. (3.14) 



ab,N,- 



Let us now calculate the matrix elements of the hamiltonian associated to 
^it- note that the only matrix element different from zero is 



{{ab,g,s)-v\HH(^'b',n,s',A)-v') 
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-cA* 

= (2nf5^^\ma,{v-v'))EXn'^Xn'iO)^^^^ (3.15) 

where 

{{ab,g,sy,vKT^r'himO) =: f^^J,,, . (3.16) 

In the calculations above we have not made any explicit assumption about the 
relative size of and L^. We are mainly interested in very heavy quark- antiquark 
bound states where small momenta can be considered as corrections, at least for the 
lower energy levels. Clearly, these bound states should be singlets since the octet 
potential is repulsive. In fact, at the level we are working, the octet states are not going 
to give contributions to the physical observables so we will neglect them in the following. 
Hence from now on colour singlets are understood and colour indeces dropped. We 
also remark that we are always working in the CM frame, even though sometimes we 
keep -y 7^ in some intermediate steps for convenience. Following standard Quantum 
Mechanics perturbation theory |]l^ we can obtain the corrected bound state energy^ 



and wave functions (states) for the lower levels. They read 

SEab,n = ""'^^"''"V aM$aM(Q)r (3.17) 

I {ab, n,s);v)F = Z]!"^ \ {ab, n,s);v)F (3.18) 

I {ab, n,s);v)F = \{ab, n, s);v) + \ {ab, n,s);v)^^'^ + \ {ab, n, s); ■u)'-^^ + ... (3.19) 

I (a6,n, .);.;)« = ^^^^^^^!^G'^(E„,,„) J5.(0)|0) (3.20) 



\{ab,n,sy,v)^'^ = E I {ab, m, s); v) "°^^^-^-") Eab,nKbAO)^aUO) ""^''"^ 



(3.21) 

Zn^l + ^ '^"ri^°''"^ l-^'^M'^aM(0)P (3.22) 
where both continuum and bound states are included in the sum in ( |3.21 ), (|3.1 



denotes the physical normalized state (with low momentum corrections) and 

^'(^) ■= — mr-- (3-23) 

z — + le 



^The correction to the bound state energy was found to be zero in [g[ because the existence of L^^ 
was not known. 
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; I dSe*^"^(o|r { 7,^.(0)} |o) =: n,,(E,,,„)tr(r'r) (3.24) 



Pn = imab,n, 0) , mab,n ■= rUab + Eab,n ■ 

We should stress that in the last two equations there is only small momentum dynamics. 
High energies may come from the external bound state energy insertion. 

Some comments are in order. Notice first that for / 7^ (angular momentum) the 
wave function (state) and the energy remain unchanged. Notice also that the s-wave 
state does not receive contributions from / 7^ states either. The previous statement 
is true due to the fact that the momentum wave function at zero momentum for / 7^ 
is zero. This means that the new interaction does not couple / = states with / 7^ 
states. This result would change if we kept further terms in the effective lagrangian 
(see ( |2.16| )) but, of course, these contributions would be sub leading. 

Let us next calculate the decay constant. In order to do it we split the current as 
in the last section. The soft current only gives a contribution with the low momentum 
states g in the same way as the hard current only gives a contribution with the modified 
Coulomb bound states. Thus, we obtain 

{0\4%0)\{ab,g,sy,v) = (0| jS(0)|(a6, s); = J-^^^^U^^ (3.25) 

{0\4''miab,n,sy,v) = (0|jfr(0)|(a6,n,s);t;) = -tr (FT) y^$«,,„(0) . (3.26) 
Finally the decay constant reads (changing to relativistic normalization) 
{0\4''iO)\iab,n,sy,v}p = -tr (FT) vWjV^$aM(0) 

1 AN, dEat,n ' ' ' ^ ^' 2iV, ' $afe,n(0) 



1 + E <^a^r.(0)l>:,,„(0)— I } . (3.27) 

mj^n ^ab,n ^r, 



Finally let us obtain the bilinear currents at zero recoil. For that we need to know 

J^^(O) = g'Fg^(O) , (3.28) 

{{ac,N\s')-v\J^\{))\{ab,N,s)-v) . (3.29) 

In order to deal with them we need to perform the splitting between large and small 
momentum. However this current cannot be in general splitted in two terms. We have 
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mixing between large and small momentum. Fortunately, the mixing terms disappear 
if both inicial and final states have the same velocity. This will not longer be true for 
non-zero recoil matrix elements. Thus, we obtain 

{{ac,n\s'yMJTm{ah,n,s)-v) = Jl^llTi f (3.30) 



{{ac,g',s'y,v\4'^miab,g,sy,v) =: - "^^J^^^ fZ, (3-31) 
where we have used 

)a2a4,{^ )aia3 ~ 2(p_|_)q,203 (P-)oi04 ■ (3.32) 

s 

The remaining possible matrix elements are zero. Notice that fab% = ^g'g because of 
baryonic charge conservation. 

The physical matrix element reads (again with relativistic normalization) 

F{{ac,n,s');v\J^''{0)\{ab,n,s);tf}F = -^/mab,nmac,n''tr{T''T'T) 
{ r d?k ~ ^ ~ 

^ liacAEac^'^ Eab,n) E,,^^, E,b,nKcA^)^^UO) (3-33) 



Eah,n — Eab,m J IJ- (2vr) 



Eac,n' ~ Eac,m (271") 

where 

j ti^xirfV2e^^n^^-^ie-^^"-^'^(0|T{j,';^.,(xi)j;^r(0)J"^ |0) 

=: I{acAEac,n',Eab,n)tr{r'T'T) (3.34) 

We can easily check that orthonormality is fulfilled when 6 = cQ. We expect the last 
statement to be true since spin symmetry relates the matrix element with the baryonic 
charge when b = c. 



^This was not always the case for the resuh given in 
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F{{ab, n', s');v\4\0)\{ab, n, s)-v)f = -mafe,„tr(r'rT)5„,„, . (3.35) 

Before finishing this section let us make some remarks. Both correlators (|3.24|) and 
(|3.34| ) should be small quantities for perturbation theory to hold. This is the case if fi 
is small against the typical momentum in the Coulomb interaction (i.e. fiaab,n << 1, 
where aab,n = n/fiabO: is the Bohr radius). It constraints the possible applications to 
the lower energy levels. On the other hand Aqcd << /i should hold so that the low 
momentum dynamics is not strongly affected by the cut-off. 



4 Cut-off independence 

Our results in the last section may look like strongly cut-off dependent. We have 
two sources of cut-off dependence. On the one hand we have a cut-off separating small 
momentum gluons from large momentum gluons. This cut-off is the responsible for 
the absence of Coulomb interation in L^. It has been mentioned at several instances 
but it has never been written down explicitly in the formulas. This cut-off dependence 
has been analysed before [|l^ so we shall ignore it in the following. On the other hand 
we have the cut-off separating large and small relative momenta. It plays the role of 
an infrared cut-off in the perturbative Coulomb wave function (large momentum) and 
the role of an ultraviolet cut-off for the small momentum contributions. We prove in 
this section the cut-off independence to the desired order (/x^, Aq^j^) of this last cut-off. 
This is crucial to ensure that our approach respects colour SU (3) gauge symmetry. It is 
important to use the same cut-off procedure in both large and small momentum regions 
in order to neatly cancel the cut-off dependence. We use a hard three momentum cut-off 
for convenience, as we have done in the previous sections. 

First of all, let us study the cut-off dependence in the low momentum correlators we 
found in the last section. Although they are non-perturbative objects we can always 
perform a perturbative calculation in order to see how they depend on the cut-off. 

Let us start by (|2.13|) ( which is formally equal to NRQCD). For ( |3.24|) we obtain 
at the lowest order (in the CM frame, v = 0) 



1 , /I +X 

— In 



2x Vl — X 



x= , ^ (4.1) 
2/i„,(fco + ^e) 



Let us consider two limits. 
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In the limit x >> 1 (i.e. near threshold) it reduces to 

(4.2) 



n.,(fc°) ^ 



1 . 



and no pole appears. In the limit x << 1 ( |4.1| ) reduces to 



This expression is going to be important in the following. We stress that (O) is fiab 



independent, and amounts to drop the 1/m terms in (|2.13|) which is nothing but the 
HQET for quarks and antiquarks. Let us now look for the physical situation we are 
interested in. Thus, we take = E^.„ and we obtain Ixl = — — , but this is 
nothing but the parameter we need to keep small so that the small relative momentum 
contributions are subleading, and hence our expansion makes sense. In the following, 
we always consider that we are in the limit |a;| << 1. 
For (p.34|) we obtain {x « 1) 

At this point we would like to stress that in the limit x << 1 the same perturbative 
results are found using HQET. This is going to be determinant in the next section. 



and (|4.4| ) make explicit the UV cut-off dependence coming from the small 
relative momentum region. Let us next go on with the IR cut-off dependence coming 
from the large relative momentum region. 

Let us then study the cut-off dependence of the wave function (for simplicity we omit 
the flavour indices). In order to do it we solve the wave equation ( |3.11| ) perturbatively 
in /i. n labels continuum or discrete spectrum. Because of the radial symmetry, we 
can write (we follow HTSl) 



^n,l,m{P] yU) = Fnj{p] fi)Yi^m{p) (4.5) 

where Fn^iip] jj) fulfils 

Fn,i{p] fi) / qdqFn,i{q; f^)Qi{— ) = En{n)Fn,i{p] fi) p > fi (4.6) 

Trp Jfi 2pq 
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Qi{z) = \f dx^ (4.7) 

Z J-l z — X 

and Pi is the Legendre function of the first kind. We stress that we are interested in 
Fn,i (p; A*) ioT p > fi only, although in the intermediate steps it is going to be defined 
over all p > values. Now we perform a cut-off parameter expansion and we work as 
in usual quantum mechanics perturbation theory where we demand the corrections to 
be orthogonal to the leading result. 

oo r oo r 

r=0 ' ■ r-=0 

We also expand the cut-off in the integral. 

/ qdqF:,{q)Qi{^^^) ^ h^.ip, f^) ^f^hZip)^ • (4.9) 

On general grounds we can see that the corrections to the Coulomb wave function 
and energy go hke 0(//^'+^), therefore, as expected, we can neglect the I ^ states 
since their contributions are subleading. 

At leading order we obtain the standard Schrodinger equation with a Coulomb 
potential with no dependence. Furthermore, for the following terms in perturbation 
theory we obtain 

El^El^^, |.i(p-)=|.2(p-)^0. (4.10) 
Finally to third order we obtain 

B» = -£„M, (4.11) 

(4.12) 

We have not yet normalized the cut-off dependent wave function, as we can see from 



Therefore, we must change 



$„(p;/.) ^ ^MlA (l + l*n(0)P^) . (4.14) 
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([4.11| )-( p:.14| ) provide the explicit IR cut-off dependence from tlie large relative mo- 



mentum region. 

We have obtained the explicit cut-off dependence to the desired order /i^ in both 
large and small momentum regions. Now we will see they match properly, that is, the 
observables are cut-off independent. In fact what we will see is that the physical states 
(|3.18| ) themselves are already cut-off independent. In this way we prove the cut-off 
independence for any observable. 

Consider first the bound state energy 

E^ = E:^ + 6E^. (4.15) 



The cut-off dependence of En is given by (^.8|), ( [4.10[ ) and ( [4.11| ), whereas the cut-off 



dependence of 5En is given by (|3.171 ) and (|4.3|) . One can then easily check that E^ is 
cut-off independent. 



Consider next the state \{ah,n, s)\v) p in (|3.19|) . Recall that the first and last term 



on the rhs belong to the large relative momentum region whereas the term in the 
middle belongs to small relative momentum region. Let us keep apart for a moment 
the explicit cut-off separating these two regions in the relative momentum integrals. 
The remaining cut-off dependences of the first term are given by ( |4.5|) , (|4.8| ) , ( [4.10|) and 



(|4.12| ), while for the last term are given by ( p.21j ) and (^4.3| ), which cancel each other. 



It remains the UV cut-off dependence coming from ( |3.20| ) (which has been already 
studied in 0) and the explicit IR cut-off dependence coming from the integral over 
relative momentum in the first term of ( p.l9| ) (see (|3.1|) ), which we kept apart for 
a while. Recall that the wave function in the first term of ( |3.19| ) is, except for the 



normalization factor ([4.14|) , the Coulomb wave function since we have already cancelled 



the cut-off dependences coming from (|4.12| ). Let us next calculate (|3.20|) perturbatively 
at lowest order. It reads 

\{ah,n,sW^ = f ^t."(Pi)r.t;^(P2)4bi)4UP2)|0) 

^aUk)u''iPi)^sV^{p2K{Pi)d%{p2m- (4.16) 



The second equality holds at the order we are working at. Notice finally that this is 
nothing but the piece we need to add to the first term of ( p.l9| ) in order to obtain a 
relative momentum integral independent of the cut-off. Finally, the cut-off dependence 
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of the normalization in ( |4.14| ) and of ( |3.22| ) also cancel each other in ( |3.18| ) (again 
taking into account ([4.3|)). 



We have thus seen that at the level of physical states we are able to prove the 
cut-off independence. The cut-off independence can also be checked explicitely in the 
observables (|3.17|) , (|3.27|) and (|3.33| ). This demonstrates that the HQET ultraviolet 



behavior cancels the NRQCD infrared behavior in Coulomb type bound states, which 
garanties that we have performed a proper matching between large and small relative 



momentum. This issue has also been pursued in p, 15 . 



5 Evaluation of the low momentum correlators 

In section 3 we learnt how to parametrize the possible non-perturbative contribu- 
tions in the small relative momentum region in terms of two low momentum correlators 
( (|3.24|) and (|3.34| )) with external Coulomb bound state energy insertions. It is remark- 



able that these contributions only exist for s-states. At the beginning of section 4 we 
also saw that the kinetic term, which is suppressed by a mass invers power, can be 
safely neglected in the correlators we are interested in, and hence we can use HQET 
for quarks and antiquarks to discuss their properties. 

The HQET for quarks and antiquarks enjoys a U{4Nhf) symmetry which breaks 
spontaneously down to the U{2Nhf) ® U{2Nhf) Isgur-Wise symmetry |1T2| . 

Let us first analyse the consequences of the unbroken U{2Nhf)®U{2Nhf) symmetry. 
In fact the spin symmetry which is included in it has already been used in (p.24|) and 
(|3.34| ). The flavour symmetry implies the following 



fab,g — fg , fac%b — ^ ■ (5-1) 

Therefore we get 

^ab{Eab,n) = ^l{Eab,n) , ^ac,ab{Eac,n' , Eab,n) = ^2{Eac,n' , Eab,n) ■ (5.2) 



The correlators (|]2|) and (|3^ ) are thus given in terms of two unknown universal 



(flavour independent) functions Hi and H2. But if we go further, using flavour number 
conservation together with flavour symmetry, we obtain /^'^ = 6gig for any flavour. 
From that it follows that if Hi is known for any energy insertion we can obtain H2. 
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Explicitely they read 

Ili{Eab,n) = ' ^2{Eac,n', Eab,n) = 



n TP 7? ' u,u,ii/ ^ ^ O fP WW W 

^ ^ah,n g ^ ^ac,n' ^g ^ah,n 

(5.3) 

These low momentum correlators can be further specified at least in two situations. 

i) Eab,n » ^QCD , {mq -> oo , « Small) , 

ii) Eab,n « J^QCD , (ruQ large , a ^ 0) . 

Notice that situation ii) is conceivable if a is very small since so far we have only 
assumed that the invers Bohr radius is much bigger than Aqcd and the energy is 
suppressed by a factor a with respect to the formeiQ ^ 

In the situation i) the operator product expansion holds. If we carry it out for the 
low momentum correlators we just obtain ( [4.3|) and (|4.4|) . Their cut-off dependence just 
cancels the cut-off dependence from the large relative momentum region, as we saw in 
section 4. Hence, we conclude that there are no new non-perturbative contributions in 
this situation, thus confirming the fact that the VL contributions from the condensate 
are indeed the leading non-perturbative effects in the mq —>■ oo limit]^. This result 
follows from the observation that there is no local gauge invariant object that can be 
built out of Dq alone. We have explicitely checked it for lower order terms. 

In the situation ii) we are in the low energy regime of the HQET. In this regime 
it is important that the HQET with quarks and antiquarks with the same velocity 
undergoes a spontaneous symmetry breaking of a U{ANhf) symmetry down to the 
Isgur-Wise symmetry U{2Nhf) U{2Nhf), since the Goldstone modes associated to 
the broken generators dominate the dynamics. The Heavy Quark Hadronic effective 
lagrangian describing the Goldstone modes was worked out in , where the correlators 
( |3.24| ) and ( |3.34| ) were also calculated. Using those results we obtain 



^In practise we must remember that a should better be substituted by the running coupUng 
constant at the quarkonium scale, which is in fact an implicit function of mg and Aqcd- 

^In , the bound state energy Eab.n was understood as giving rise to a residual mass for the heavy 
quark and antiquark in the Heavy Quark Effective lagrangian, which was latter on subtracted. That 
definitively obscures its actual role, which eventually led to some confusion: in |^ the situation ii) 
was not allowed whereas the Heavy Quark Hadronic lagrangian was used for situation i) , which is not 
correct. 

^This point was not properly specified in p]. 



21 



IlacAk", k') = Y (pT^ (^^fT^ ' (5-5) 



+ ^ (5-6) 



2 2 67r2 

where ffj is cut-off independent. Notice that in this situation all non-perturbative 
effects in the small relative momentum region are parametrized by a single nonpertur- 
bative constant which is spin and flavour independent^. This is a non-trivial conse- 
quence of the U{4Nhf) symmetry. The fact that the latter is spontaneously broken 
down to U{2Nhf) ® U{2Nhf) allows us to know the Green function behaviour at low 
energy insertion with a single nonperturbative constant since no mass term appears in 
the pole. All the spin and flavour dependence is explicitely known in the observables. 

However, caution must be taken in the situation ii). This is due to the fact that, 
in this situation the standard evaluation of non-perturbative contributions in the large 
relative three momentum region coming from ( |2.15| ) becomes unreliable. Let us briefly 
recall the two approximations involved, namely the multipole expansion and the adi- 
abatic approximation. The first one is an expansion in Aqcd over the invers Bohr 
radious, which has also been assumed to hold throughout. The second one requires 
the time evolution of the soft gluon fields to be slow in comparison with the energies 
involved in the Coulomb spectrum. This requirement is in fact the opposite of situation 
ii). Thus we are in the unfortunate even that when we have an excelent parametriza- 
tion of the non-perturbative effects in the small relative momentum region ( (^.4|) and 
(|5.5| )) we loose control of them in the large relative momentum region. 

Nevertheless, we envisage a situation where the parametrization ( p. 41 ) and ( |5.5| ) may 
be useful. Recall that although the parameter controlling the adiabatic approximation 
and the parameter controlling the expansion in the hadronic effective lagrangian are 
both of order Aqcd, they need not be exactly the same. The former was shown to be 
< DFDF > / < FF > in [0 and let 27r/| be the latter. Suppose then that 



Eab,n > 



< DFDF > 
<FF> 
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Notice also that although at first sight the contributions obtained by substituting (x4) and (5.5) 



in (3.17), (3.27) and (3.33) look like more important than those from the condensate when mg — > oo, 
they are actually not so since the smallness of a required in situation ii) mantains the condensate 
contribution dominant. Some statements made in H implying the opposite must be corrected. 



22 



^aM<2vr/|. (5.7) 

In such a situation it would be reasonable to use both the adiabatic approximation in 
the large relative momentum region and the hadronic effective lagrangian in the small 
relative momentum region. Some bottomonium, charmonium, and presumably 
states may well be considered in the situation ( |5.7| ). However, the mass of the b quark 
and mainly the mass of the c quark are not large enough to allow for a straightforward 
application of our formalism to phenomenology. Relativistic and radiative corrections 
are in general important and this is also so for the non-perturbative corrections due to 
the gluon condensate All them must be taken into account. 

Let us next discuss the expected size of our contributions. It is not our aim to 
present a full-fletched phenomenological analysis in order to extract from the data, 
which would definitely be premature as it should be clear in the following discussion, 
but just give reasonable estimates of the expected magnitude of its contributions. For 
simplicity, we will concentrate on the mass corrections. 

We start with the bottomonium system where our formalism is expected to apply 
for the lowest lying states |TU|. We proceed as follows. First of all, we fix and 
using the experimental data and the available theoretical results while ignoring 
the contribution 6Eab,n in (3.17). Then we estimate the size of 5Eah,n by letting ffj 
run within values of the order of Aqcd- We should keep in mind that although we will 
take positive for definiteness it can also be negative. We extract nib and a^,^-^ from 
the selfconsistency equation abb^iictsia'j^bi)) — ^bb,i , and the T(ls) mass. We use the 
following equation to fit the latter 



where 



mr(is) = 2m;, + A2 + A3 + Avl 
-1 mbCfa{a^^\) 



A2 = -2mb 



A, 



-2mi 



CjPoa (afefe_i)a(afefe^ 



Stt 

Avl = rrib- 



In 



+ 1-1E 



(5.8) 
(5.9) 



(5.10) 
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We have taken the formulas above, which include relativistic, radiative and the VL 
non-perturvative corrections, from |^. We allow for different values of Aqcd and 
give the relative weight of each contribution in the table. 

Let us next assume that we are in the situation (5.7). As mentioned before, this 
may well be the case for the T(IS'), T(25'), Xh(lP), J/ip (and rjc) and (and B*). If 
we let ffj^ run between the values 

Eat,n < 271 < a-l^ (5.11) 

we can give an estimate of 5Eab,n- If we allow fij between 100 — 150Mel^, our results 
turn out to be quite stable under values of Aq{t^s between 200 — 300 MeV. We obtain 

- 9MeV < 5Ehb,i < -2MeV (5.12) 
where the explicit expresion used for calculating 5Eab,n reads 

5E,bn = -4/ia, ^^4^ i (5-13) 

A'cC/a,(a„^„) \2nab) 

Although the smallness of the result above is discouraging at first sight, it justifies the 
procedure used and makes it self consistent. 

For n = 2 we obtain 



- bbMeV < 6Ebb,2 < -l5MeV (5.14) 

Recall that only the s-wave states receive this correction. If the sign of ffj was negative, 
the signs above would be reversed. This would help to understand the mass difference 
between the Xb(l-P) and the T{2S). 

Let us finally give some estimates for 6Ecc,i and SEbc,i corresponding to the J/\E' 
(and rjc) and the Be (and B*) ground states. We have taken the mass of the charm 
rric = 1570Me\^ as given in reference 0. For the J/\I' we find 

42MeV 

taking Aq^^d = SOOMel/, a'^^ = 8A8MeV and fj^ = 150MeV. For the B^ we find 

5Ebc,i ~ -23MeV 

^^However, we have not taken into account the contributions of order 0{a'^,a^) given in since 
the complete calculation at this order is still lacking. 
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taking Kq^cd = SOOMeV, a^^\ = lOlSMeV" and f]/^ = 150MeV. 

The above contributions for the energy shifts are, on the one hand, small enough 
to make us confident that our results are under control and, on the other hand, large 
enough to hope for its eventual observation. However, it is important to realize that 
the VL contributions are excedingly big for T{2S), Xb{^P), J/'ip (and rjc) and Be (and 
B*). We suspect that the framework used so far to calculate the VL contributions in 
the large relative momentum region is not appropiated for these states. We believe 
that in order to make realistic QCD-based predictions for these states one should 
devise a reliable approximation in the large relative momentum region to deal with the 
situation (ii) above , namely invers Born radius and energy larger and smaller than 



Aqcd respectively. Work in this direction is in progress |2Q 



6 Conclusions 

We are confident that the theoretical framework above is going to be useful for an 
eventual QCD-based formalism attempting to encompase situations where the Coulomb 
energy is large (small n) and situations where it is small (large n) with respect to Aqcd 
in heavy quarkonium. Even more, this formalism could also be useful in order to obtain 
explicitely the perturbative Coulomb corrections to the non-perturbative heavy quarks 
bound states (large n). 

Our formalism is clearly inspired by the Wilson renormalization group approach. 
We separate the fields into large and small momentum components by an explicit cut- 
off, and work out what the effective action for the latter is. However, there is an 
important point which makes our formalism rather peculiar: integrating out the large 
momentum components does not give rise to local counterterms only. There is non- 
trivial physics in the ultraviolet, namely Coulomb type bound states. As far as we 
know, this is the first example of a Wilsonian approach where effects due to bound 
states have been taken into account. 

Let us finally summarize the main contributions of this paper. Elaborating on the 
ideas first presented in , we have produced a detailed derivation of the effective the- 
ory governing the small relative momentum degrees of freedom in heavy quarkonium. 
In particular this includes an interaction term, which had been overlooked before, that 
leads to a few corrections in the observables. We have proven the cut-off independence 
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of the formalism. We have also discussed in detail when non-perturbative contribu- 
tions which cannot be expressed in terms of local condensates arise, namely when a 
description in terms of a Heavy Quark Hadronic Theory is adequate. Our prelimi- 
nary estimations suggests that these contributions lead to energy shifts of a few tens 
of MeV. Unfortunately, more theoretical work is necessary to establish them from the 
data. This is mainly due to the lack of control on the non-perturbative effects in the 
large relative momentum region of most of the systems where our approach should 
apply, namely T(25'), Xbi^P), J/i^ (and rjc) and Be (and B*). Work in this direction 
is in progress pO| . 
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A2{MeV) 


A^{MeV) 


AyL{MeV) 


mb{MeV) 




200 


-314 


49 


25 


4850 


1234 


250 


-376 


61 


18 


4879 


1354 


300 


-440 


74 


13 


4906 


1468 



Table 1: We display A2, A3 and Avl defined in (5.10). The last two columns give our 
results for rrib and a^\. 



